The derivative expansion of the renormalization group 
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By writing the flow equations for the continuum Legendre effective action (a.k.a. Helmholtz free energy) with 
respect to a particular form of smooth cutoff, and performing a derivative expansion up to some maximum order, 
a set of differential equations are obtained which at FPs (Fixed Points) reduce to non-linear eigenvalue equations 
for the anomalous scaling dimension r\. Illustrating this by expanding (single component) scalar field theory, in 
two, three and four dimensions, up to second order in derivatives, we show that the method is a powerful and 
robust means of discovering and quantifying non-perturbative continuum limits (continuous phase transitions). 



In this short review I want to stress what I 
believe to be the essential steps involved in the 
scheme briefly described above[2,3], and compare 
and contrast other analytic approximation ap- 
proaches to the RG (Renormalization Group). 
Particularly in the two dimensional example [4], 
the method is seen to be in several ways more 
powerful in fact than any other approximate ap- 
proach, analytic or otherwise. This is also an ex- 
cellent opportunity to refer to some papers which 
I was unaware of at the time of writing refs[l-3]. 

The first point to make is that any continuum 
approximation method, if it is to be of general ap- 
plicability, must be formulated within the frame- 
work of Wilson's RG[5]. This is true, even if one 
looks only at perturbation theory to all orders [1], 
but is particularly obvious when the theory at 
the UV cutoff scale Ao, where the bare action 
Sa is formulated, is also strongly interacting. In 
this case the usual power counting arguments, to 
determine which are the renormalizable (i.e. rel- 
evant) interactions, breaks down and the defini- 
tion of the Lagrangian itself becomes part of the 
calculational problem. Breaking this circularity 
was one of the many successes of RG viewpoint. 
In terms of this one should think of the space of 
all possible local interactions. The theory is then 
renormalizable if it is describable in terms of self- 
similar evolution in some subspace of parameters 
under the flow of some low energy scale A « Aq . 



In the way this is usually formulated one thinks 
of an effective action T\ designed to work with 
a UV cutoff A in such a way that it reproduces 
exactly the same physics as one obtains with Sa 
(and one often thinks of rescaling the cutoff back 
to its original size after the "blocking" ) but these 
concepts are not really necessary. 

The next question to ask is: What are the pos- 
sible approximations in this approach? The ob- 
vious approximation is to truncate the space of 
interactions to just a few operators. This direc- 
tion has been extensively pursued over the years, 
especially within real space RG[6], and the results 
can sometimes be very accurate - and sometimes 
not. Although by no means as sophisticated yet, 
approximations by truncations of the continuum 
RG also have a long history [5, 7]. However, if the 
field is strongly interacting, the couplings of all 
the interactions <fp ,n will be more or less the same 
size, so these methods need a lot of care. (Reasons 
for some limited accuracy are given in ref.[3]). To 
have a chance of finding a robust method of gen- 
eral applicability we must keep at least infinitely 
many such operators. 

We will do so by using a <9E (Derivative Expan- 
sion) up to some maximum order. (Maybe there 
are other ways). Thus to 0(d 2 ) for example we 
write 



Jd D xV(<P)+ l -{d^YK(<P) 



(1) 
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with general functions V and K, and ignore 0(d 4 ) 



terms and higher. This expansion seems partic- 
ularly natural, and one might hope that it has 
general applicability in the sense that if the higher 
derivative terms are not small, then the descrip- 
tion in terms of a single scalar field is probably 
itself inappropriate and indicates that other de- 
grees of freedom should be introduced. 

So far we have dealt with rather unsubtle re- 
quirements, but there is another requirement, 
which determines the form of the cutoff, and to 
highlight this it will help to consider the simplest 
sort of self-similar evolution, namely when Ta 
is actually invariant under lowering A, i.e. FPs 
(Fixed Points) of the RG. Since these occur at 
zeroes of j3 functions, we generally expect only a 
few such FPs, and that at each FP all parameters 
will be dynamically determined. Independence of 
A implies scale invariance, and thus a massless 
theory. As is well known, this scale invariance is 
realised only if one assigns the field the dimen- 
sion [4>] = dj, = 7}(D — 2 + rj), where D is the 
space-time dimension and r\ is the anomalous di- 
mension. A must be introduced in a way which is 
equivalent to an IR cutoff C(p, A) [4]. Introducing 
it into the partition function as 

Z[J] = JV4> exp{-^.C- 1 .0-5 Ao [0] + J.0}(2) 

it is then easy to write down a flow equation: 
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Z. Now we rewrite all 



2 8.J- dA -8J d 

quantities in terms of dimensionless quantities, 
using A, so that everything will be independent 
of A at a FP. (From (2) by dimensions, we require 
[C] = 7j — 2). It is also very helpful to rewrite it 
as an equation for the Legendre effective action 
T [</>]. There are a number of reasons for this[l, 
2,4]; I will mention one later. Having done this 
we obtain an equation for T[<j>] with the following 
features [2]: It has 5 2 TI54> 2 (as a consequence of 
the two J derivatives), and n appears explicitly 
now as a parameter in the equations. 

Consider what happens when we look for FPs 
at lowest order in the <9E, O(d ), where we just 
set K = 1 in (1). Since this corresponds to throw- 
ing away all momentum dependent corrections, it 
is clear that n = in this approximation. Substi- 
tuting this into our flow equation, it follows from 
the previous description that one obtains a (non- 



linear) second order (ordinary) differential equa- 
tion for V((f>). Such an equation has generically a 
two parameter set of solutions. If these all looked 
acceptable it would indicate that we would need 
further information and cannot just examine the 
FPs directly in this way. Fortunately some magic 
occurs. It turns out that (generally) very nearly 
all solutions are unacceptable because they are 
singular at some finite real <j> = <j> c [2-4] . The few 
non-singular solutions are approximations to the 
physical FPs. Effectively O(d ) equations have 
been considered before[8,9] with similar conclu- 
sions^, 10-12]. If we fix V'(0) = 0, considering 
only (j> "B- — <f> invariant theories, then the value 
of <j) c can be plotted against, say, V(0). In this 
way we can scan through the infinite dimensional 
space of all possible potentials V ((/>), looking for 
continuum limits (<f) c = oo) [3]. Clearly this is 
much more than is possible with other approxi- 
mate methods. However, only the expected FPs 
have been found at this level. In particular only 
the Gaussian FP in 4 dimensions[ll,3], this and 
the Wilson FPs in 3 dimensions[10-12,2,3], and in 
2 dimensions only critical sine-Gordon models [4]. 
There are very strong conjectures for the exis- 
tence also of an infinite set of multicritical FPs[13] 
in 2 dimensions, but this requires n > 0. I will 
come back to these, but first I want to show that 
there is another way of seeing that there are only 
a few FP solutions: On the one hand V'(0) = 0, 
and on the other hand, if the cutoff A drops out 
we have, in the limit <j> — > oo, V(4>) oc <f> D / d *, 
by dimensions. This gives us now the needed two 
boundary conditions which indeed are satisfied by 
the FP solutions (providing d^ ^ 0). 

It is useful to extend this argument to 0(d 2 ) 
and higher. Thus consider the FP equations one 
obtains now by using (1) with K a fully fledged 
function. Just as before, the structure of the flow 
equations means that we get second order ordi- 
nary differential equations, now a coupled pair of 
them. Now we expect two extra parameters, but 
we also have now two more constraints, coming 
from symmetry and dimensional arguments on K . 
In this way, to all orders of the <9E, the coefficient 
functions are determined for each FP. However we 
have been ignoring n, which can now be non-zero, 
and if the method is to work at all it must also get 



determined somehow. It gets determined in the 
exact (= 0(9°°)) equations via a reparametriza- 
tion invariance[14] of the flow equations, involv- 
ing a rescaling of the field[14-16]. The point is 
that this extra invariance means that we can fix 
an extra condition, which now overconstrains the 
equations, turning them into (non-linear) eigen- 
value equations for r\. (E.g. use it to normalise 
the kinetic term as A'(0) = 1, as in fact was done 
implicitly at 0(d )). The problem is that this 
invariance is generally broken by the approxima- 
tion scheme[15], including the <3E[17], and one is 
reduced to heuristics for trying to find a 'best' 
value for ij[15— 17]. Actually there is a, probably 
unique, form of cutoff for which the <9E preserves 
this invariance: C(q 2 ) oc q 2K , for some integer 
K > D/2 - 1 [2]. This only cuts off the Legen- 
dre equations, and only if this inequality is satis- 
fied. It works because the invariance is then just 
another scaling symmetry with a (non-physical) 
set of scaling dimensions. Clearly the <9E will 
converge best for the smallest value of re[l,2]. 

For the operator spectrum, nothing new is 
needed. Linearising the flow equation around the 
FP solutions, we have again precisely the right 
number of boundary constraints, from symme- 
try, dimensional analysis and the normalisation 
allowed by linearity [2]. 

It remains to state that the results of this 
method are so far very encouraging. In particular 
to 0(d 2 ) one still finds only Gaussian and Wilson 
FPs in three dimensions, but now with estimates 
for critical exponents which are within six times 
the error of the best present estimates [2]. In two 
dimensions, a constrained search finds only the 
conjectured r\ > multicritical points [4]. For the 
first ten such scalar field theories (critical, tricrit- 
ical, quadricritical, ..., undecimcritical) the expo- 
nents and dimensions of up to the first ten op- 
erators agree within 33% with expectations[13] 
over two orders of magnitude, with a weak im- 
provement with increasing multicriticality and di- 
mension. (We may also have found "shadow op- 
erators"^]). These facts are particularly signif- 
icant, since all other present approximate non- 
perturbative methods, lattice Monte Carlo, re- 
summations of weak or strong coupling perturba- 
tion theory and the e expansion, become rapidly 



impracticable with increasing multicriticality or 
dimension [4]. 
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The two dimensional undecimcritical point. 
V(4>) multiplied by 33 for display purposes. 
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